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Programa 


1. Introdução à Teoria Geral dos Sistemas. Noção de sistema. 
Principais características dos sistemas. Hierarquia de 
sistemas. Sistemas dinâmicos. 

2. Modelização de sistemas por equações diferenciais. 

3. Função de transferência e dinâmica de sistemas. Resposta 
temporal de sistemas dinâmicos 

4. Representação no espaço de estados. Estados singulares e 
linearização. 

5. Dinâmica de sistemas não lineares. Sistemas caóticos e 
fractais e sua aplicação à Biologia e Medicina. 


Simulação computacional de casos de estudo em ambiente 
Matlab/Simulink. 
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Modo de funcionamento 


Aulas teóricas de discussão da matéria (conveniente 
estudo prévio, materiais de estudo já disponíveis) 


Aulas práticas de implementação computacional e 
simulação em Matlab/Simulink 


Modo presencial 
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Capítulo 1 


Noções Fundamentais de 
Sistemas 
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1.1 Sistemas, elementos, relações 


1.2. Retroação ou realimentação (feedback) 

1.3 Ambiente, sistemas fechados, sistemas abertos e homeostasia 
1.4 Entropia e neguentropia 

1.5 Estado e espaço de estados 

1.6 Autopoiese, adaptação, cibernética 


1.7. Uma hierarquia de sistemas 


1. 8. Sistemas dinâmicos: noção formal 


1.9. Conclusão 
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1.1. Noção de sistema 


Exemplo 1 
( 
= 3 esa 
Vaso com água Lamparina Fósforos 


... Separados não constituem um sistema. 
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Exemplo 1 


Elementos: o vaso, a água, a lamparina. 


Comportamento: aquece a água 


Atributo: temperatura, 2C. 


- Separados, os elementos não constituem um sistema. 
- montados, constituem um sistema. 


- a montagem consiste no estabelecimento de relações 
entre os elementos (proximidade, transmissão de 
calor, etc.). 


- efeito de sinergia: o todo é maior do que a soma das partes. 
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Exemplo 2 


Um automóvel Elementos: as suas numerosas partes. 


Comportamento: posição, velocidade, 
aceleração. 


- desmontadas: não constituem um sistema. 
- montadas: constituem um sistema. 


- à montagem consiste no estabelecimento de relações 
entre os elementos. 


- efeito de sinergia: o todo é maior do que a soma das partes. 
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Exemplo 3 Um computador 


Comportamento: velocidade, capacidade 
de armazenamento, temperatura, ... 


Elementos: as suas numerosas partes 
(CPU, discos, teclado, monitor, etc.). 


desmontadas: não constituem um sistema. 


montadas: constituem um sistema. 


a montagem consiste no estabelecimento de relações 
entre os elementos. 


- efeito de sinergia: o todo é maior do que a soma das partes. 


(OADC/DEI/FCTUC/LEB//MCPF/DSBF/2022/Capl1 10 


Exemplos 4 


Sistema de ar condicionado 


Temperatura 
desejada 


Temperatura 
medida 


Realimentação negativa (Feedback) 


Sistema de controlo, servomecanismo, regulador. 
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Exemplos 5 


Célula animal 


Cross-Section of an Animal Cell 


Cell Membrane Centrosome 


Cytoplasm 


Rough ER 


Nucleus 
Nucleolus 
Smooth ER 
Nuclear 
Membrane 
Vacuole Ribosomes 
: R Golgi Body 
ia OEnchantedLearnina.com 
( de ) 
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Exemplo 6. Sistemas fisiológicos 


Uma célula 


Inhalation 
—----——-- Exhalation 


O coração 


O sistema respiratório 


( poetas 
pressure 
Rib cage Diaphragm 


O sistema nervoso 
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Exemplo 7. Sistema humano 


Uma pessoa, a sua história, os seus sonhos e as suas utopias 
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Exemplo 8. Sistemas sociais 


Um coro 
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Lo 


Exemplo 9. Sistemas económico-sociais 


Uma empresa 


Uma cidade 


. níveis de complexidade crescente. 
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Sistema: definição formal 


> Um conjunto de elementos (componentes) 
relacionados e interagindo entre si, num todo 
organizado, de modo a desempenharem 
alguma função. 


Elemento : algum fenómeno da vida natural ou social, ou a 
sua representação, que tem certos atributos que podem 
variar ao longo do tempo em consequência do seu próprio 
comportamento. 
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Não é sistema Não é sistema 


É sistema 


BOUNDARY , 


INPUT 


(De Flood and Carson) 
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Existe uma relação entre os elementos El e E2 (de um 
sistema) se o comportamento de E1 é influenciado ou 
controlado pelo de E2 e se o comportamento de E2 é 
influenciado ou controlado pelo de El. 
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Um atributo de um elemento é uma propriedade característica 
do elemento (cor, tamanho, peso, temperatura, etc.). 


Um atributo de uma relação é uma qualidade característica da 
relação (intensidade, velocidade de comunicação, etc.). 


No estudo de sistemas dinâmicos a evolução dos atributos 
assume um aspecto fundamental, dado que é essa mudança 


que exprime a dinâmica. 
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1.2. Retroação ou realimentação (feedback) 


existe quando a influência de um elemento inicial se repercute 
noutros elementos e através de uma série de relações o efeito 
dessa influência inicial acaba por se repercutir sobre o elemento 
inicial. 


Exemplo: a dinâmica de um par predador- presa: 


> Um aumento da vegetação produz um aumento explosivo da 
população de um pequeno herbívoro , 


> o que possibilita que a população de um carnívoro, predador 
do herbívoro, exploda também, 


> e em consequência a população do herbívoro diminui. 


> O aumento da população de herbívoros retroagiu 
negativamente sobre si própria através das suas relações 


com a população dos carnívoros. 
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Retroação negativa 


> Para atingir objectivos estabelecidos através de parâmetros de 
controlo (quer feitos pelo homem, quer naturais) 


> Se um sistema se afasta do seu regime estacionário, ou se 


executa uma acção de controlo, ou gera-se retroacção natural 
para reverter a situação. 


No caso do par predador - presa 


E aumenta E 
População de A população de 
herbívoros carnívoros 
aumenta 


Hi |. E 
diminui 
DOOR OO OO OO 
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A retroação existe no corpo humano em muitas situações, 


de complexidade variável. 


Cerebral cortex Regulation 
device 
Q Central 
Governing value switching 
mechanisms N, vagus -—» 
Diencephalon E B-cells Insulín 
Pancreas 
islands 


A-cells 


Medulla 
Adrenals 
Cortex 


Corticosteroids 
Adrenalín 
Insulín 
4-— Feedback 
Regulation 
value 
Blood sugar Liver storage 


Glycosensíble 
receptors 
tn 
pancreas, 

C.N.s. 


Muscle energy 
consumption 


Kidney 


Blood overflow á 
(hypothetical canto valve 
Governor 
t 
(b) Lone (Bertalanffy) 
STIMULUS MESSAGE MESSAGE RESPONSE 


CONTROL 
— >» | RECEPTOR |—y EFFECTOR | ——» 
RECEPTOR [+ | appaparys | —>[ EETECTOR 


(2) FEEDBACK 
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Regulação 
do açúcar 
no sangue 


23 


Retroação positiva 


> Também existe. 


> Pode ser desejável mas pode levar a mudanças estruturais 
e até ao colapso do sistema. 
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1.3. Ambiente (do sistema), sistemas fechados e 
sistemas abertos e homeostasia 


Ambiente, 


> Com o qual o sistema partilha apenas relações de 


entrada/ saída. 
Fronteira 


> A demarcação entre o sistema 


ENVIRONMENT 
NE 


LEIA 


e o seu ambiente. 


NARROWER SYSTEM 


WIDER SYSTEM 


ENVIRONMENT 


WIDER ENVIRONMEN 
(De Flood and Carson) 
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Sistema aberto 


> Troca materiais, informação e/ou energia com o ambiente 


> 


> 


através de uma fronteira (como a identificar?) 

Mantém-se a si próprio continuamente a construir e a 
destruir componentes (fluxos de entrada e de saída) 

Nunca se encontra, enquanto estiver vivo, num estado de 
equilíbrio químico ou termodinâmico, mas mantém-se num 
regime que se pode considerar estacionário. 


Sistema fechado 


> 


> 


> 


Não existe nenhuma relação entre os elementos do sistema 
e elementos exteriores ao sistema 

Exemplo: a química física trata das reações, das suas 
velocidades, e do equilíbrio químico estabelecido dentro de 
um reator em que se juntaram os reagentes. 

As leis da termodinâmica aplicam-se apenas a sistemas 
fechados. 
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Uma relação materializa-se através de fluxos de materiais 
(massas), de informação, ou de energia entre os elementos. 
Por isso o conceito de comunicação pode ser usado 
equivalentemente ao de relação. 


Os três pilares da ciência moderna 
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Homeostasia 


> O equilíbrio dinâmico entre o sistema e o seu ambiente, 
com fluxos que entram (influxos) e fluxos que saem 
(efluxos). 


Ex. regulação térmica nos animais de sangue quente: quando o 
sangue arrefece o cérebro “liga” algum mecanismo de produção 
de calor no corpo, e a temperatura do corpo é monitorizada por 
retroação pelo cérebro de modo a que a temperatura do corpo 
se mantenha constante. 


> A constância de um elevado número de variáveis físico- 
químicas do corpo é mantida por mecanismos de 
homeostase. 


> A homeostase (ou homeostasia) actua para manter uma 
certa ordem dentro do sistema, através de mecanismos de 
retroacção comparáveis aos servomecanismos na mecânica 
e na tecnologia. 
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Homeostasia e regime estacionário 


ko o 


: lo 'e bo 
: rº | oureut * 
OUTPUT INPUT 
INPUT . co : 
e n 0) 0 “ Io o 
e o 


(De Flood and Carson) 


Instante t Instante t+dt 


As substâncias no interior da célula não são as mesmas, 
mas a estrutura da célula mantém-se. 
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1.4 Entropia e neguentropia 


> A entropia exprime a tendência que as coisas têm para 
evoluírem para um estado de grande desordem ou 
desorganização. 


,* INK 


& 


2 
« a 
> In > de 
WATER do 
e Flood and Cars! 


> A entropia é um processo que carrega consigo a lei 
“irreversível” da degradação da energia e da matéria. 


> E uma medida da probabilidade (dos estados: um 
sistema fechado tende para um estado de uma 
distribuição mais provável - a tendência para a desordem 
máxima). 
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<> 


Trabalha para a Trabalha para a 
desordem ordem 


> Um sistema aberto, como um sistema vivo, tem que 
importar energia, informação e materiais para compensar 
a tendência para uma maior desorganização. 


> Um sistema aberto importa entropia negativa 
(negaentropia, negentropy). 
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Num sistema fechado, a entropia cresce sempre, de acordo 
com a equação de (Rudolf Julius Emmanuel) Clausius 


15 Va 


http:/Awww-groups.dcs.st-and.ac.uk/—history/Mathematicians/Clausius.html 


a 


Lei de Clausius 


Num sistema aberto, pelo contrário, a variação total de entropia 
pode escrever-se como a equação de Ilya Prigogine 


http://nobelprize.org/chemistry/laureates/1977/prigogine-autobio.html$10  http:/Avww .fortunecity.com/emachines/e1 1/86/entropy.html 
dS=d S+dS 
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Auto-organização nos sistemas vivos 


- usando materiais e energia importados 
(actividades negaentrópicas ou homeostáticas) 


d,S < 0, negaentropia importada 
dS 20, degradação interna 
Sistemas saudáveis: 
ds | 2 | d;S| 
dS+dS <0 


US su 
Aumenta a ordem e a organização !!! 


No ambiente que rodeia o sistema aumenta a entropia (degradação ambiental) 
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Estrutura 


> A forma como os elementos de um sistema se podem 
relacionar uns com os outros (topografia das relações) 


> Define o contexto no qual ocorrem os fenómenos no 
sistema 


Comportamento 


> É a perceção da sua evolução decorrente de observações 
em instantes sucessivos t,, t,, ... 


> O comportamento é orientado por algum objetivo que 
pode ser avaliado em relação a uma finalidade particular 
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1.5 Variáveis de estado e espaço de estados 


(tamanhos de populações de espécies interdependentes, temperaturas dos 
corpos, volumes de água em reservatórios, etc.) 


> Os atributos elementares relevantes de um sistema 
compõem o vetor de estado x 


X 
> O estado constitui a síntese da história do x, 
sistema: ele resulta da toda a sua evolução e 
anterior 
Trajetória de estado x 


>» Quando o estado muda, os seus sucessivos valores 
descrevem uma trajetória no espaço de estados 
multidimensional. 


> Um sistema é determinístico se os valores futuros das 
suas variáveis de estado podem ser univocamente 
determinados 


>» Caso contrário o sistema é indeterminado ou aleatório. 
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1.6 Autopoiese, adaptação , cibernética 
Autopoiese 


Propriedade dos sistemas que se autoproduzem. 
Uma célula produz os seus próprios componentes que por sua 
vez a produzem a ela própria. 


Os sistemas vivos podem ser considerados autopoiéticos na 
medida em que eles se organizam para que os seus processos 
possam produzir componentes necessários para a continuação 
desses próprios processos. 


> Perturbações (ambientais) são mudanças no ambiente 
que podem fazer evoluir o sistema para um estado de 
desequilíbrio; podem ter um impacto grave no 
funcionamento do sistema 
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Adaptação 


Desenvolve-se para lidar com as mudanças ambientais e 
é necessária para a sobrevivência do sistema em tais 
circunstâncias. 

É um tipo de comportamento que visa um objectivo. 


O Darwinismo é uma teoria da adaptação. 


Se o ambiente é constante, a adaptação não é uma 
capacidade crítica. 
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Sistemas adaptativos 


Comportamento adaptativo a curto prazo 


Para lidar com impacto agudo no sistema, este deve aplicar rapidamente 
regulação e controlo para reganhar equilíbrio 


Comportamento adaptativo a longo prazo 


Se as mudanças são crónicas, é necessário para manter a integridade 
do sistema 


> Cibernética 


a ciência do controlo e comunicação nos animais e nas 
máquinas. Descreve as leis naturais que governam a 
comunicação e o controlo em situações dinâmicas. 
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1.7. Uma hierarquia de sistemas 


Metasistema (ou supersistema) 


Um sistema que se coloca por cima de outro, numa hierarquia de 
funções de monitorização e controlo 


No humana, o cérebro (com consciência) é um metasistema 


Organização hierárquica: representação lógica dos fenómenos 
como sistemas e subsistemas 
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Reducionismo sistémico 


> Aumenta o nível de resolução de análise, vendo os 
sistemas em maior detalhe 


> O que define o nível apropriado de resolução é o 
problema que está em análise, aquilo em que nos 
devemos focalizar (system-in-focus) 
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Emergência 


As organizações hierárquicas ascendentes (cada vez mais 
complexas) revelam um fenómeno importante 


(o hino dos sistemas desde há muitos anos) 


Significa que os sistemas têm propriedades emergentes 


Diferentes estruturas apresentam diferentes emergências 
Sinergia 


Palavra que exprime a emergência de propriedades 
inesperadas e interessantes 
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Metáforas em sistemas 


> metáfora da máquina 


descreve sistemas fechados com objectivos bem definidos que 
devem ser alcançados numa hierarquia rígida de controlo 


> metáfora orgânica 


descreve redes complexas de elementos e de relações com 
transformações afetadas por retroação 


> metáfora do cérebro permite emoções 


> metáfora cultural 
inclui o rede de valores, crenças e normas. 


> metáfora política 


permite redes de interesses interatuantes que as pessoas 
perseguem e podem (eventualmente) alcançar 
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Uma hierarquia de sistemas (adaptado de Bertalanfry) 


Nível 
Sistemas simbólicos 


Sistemas sócio- 
culturais 


O Homem 


Animais 


Organismos inferiores 


Sistemas abertos 


Mecanismos de 
controlo 


Descrição e exemplos 


Linguagem, lógica, matemática, ciência, morais, etc. 


Populações de organismos (incluindo humanos); 
comunidades humanas determinadas por símbolos 
(culturas). 

Simbolismo; passado e futuro, o “eu” e o mundo; 
auto-consciência resultando como consequência a 
comunicação pela linguagem. 


Aumento da importância do tráfego da informação 
(evolução dos recetores, sistemas nervosos, etc.) 
aprendizagem; início da consciência. 


Organismos do tipo planta: aumento da diferenciação 
nos sistemas (“divisão do trabalho” nos organismos); 
distinção da reprodução das funcionalidades 
individuais. 

Chama, células, tecidos, órgãos e organismos em geral 


Termostato, servo-mecanismos, mecanismos 
homeostáticos! em organismos vivos 


I homeostasia: processo de regulação pelo qual um organismo mantém o seu equilíbrio 
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1.8 Sistemas dinâmicos: noção formal 


> quanto o seu comportamento faz com que os seus 
atributos variem com o tempo. 


Muitos sistemas são atuados do exterior, através de uma 
entrada, e atuam sobre o exterior através da sua saída. 


Combustível Temperatura 
da água 


Fronteira 
com o meio 
ambiente 


A atuação pode ser através da transmissão de energia 
(térmica, mecânica, elétrica, ...) ou da transmissão de 
informação (digital, indicação de sensores, etc.). 
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Será que o comportamento num instante de um 
sistema dinâmico depende apenas da ação que se 


exerce sobre ele nesse instante ? 


Isto é, aplicando sempre a mesma entrada, obtém- 
se sempre a mesma saída ? 
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Situação A Situação B 


Temp.inicial: Temp.inicial: 
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A que se deve a diferença no comportamento ? 


O estado inicial reflecte tudo o que aconteceu antes de se 
iniciar a contagem do tempo. 


Como se chegou a 16? 


Através de aquecimentos e arrefecimentos desde -c0 até t=0 


O sistema tem memória: lembra-se, através da condição 
inicial, do que aconteceu no passado. 


A evolução futura depende do passado (através da condição 
inicial) e da entrada futura. 
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Entrada, u, Saída, Y, 


VE PO oral os qa) 


Quantidade infinita de informação necessária para 
calcular Yo... 
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Quantidade infinita de informação necessária ... 


... aproximação por uma quantidade finita 


Um e Upo Ups 


1 kk 


A 
Yen eee eee Yo Yk1 Yk 


Passado 


Futuro 


QJUSSSLA 


A saída y, é influenciada pela entrada até m instantes anteriores 
e pelas saídas até n instantes anteriores. 
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A saída y, é influenciada por m entradas anteriores e n saídas 
anteriores 


Hipótese simplificadora: f é uma função linear 


... equação linear de diferenças 


(notação equivalente) 
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Tempo discreto versus tempo contínuo 


Como exprimir a memória em tempo contínuo ? 
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TT ÃÇTÃÕÇTÃT er t=kT 


Exemplo : : : : : : 
O 123 4/56 
A discretização do eixo temporal 


y69)-yk-D=uk-n > SELO upon) 


SeT=1 > PO uk) 
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Aproximação da derivada pela diferença de 1º ordem para 
a frente: 


VD qn ss SO 


T E" =ul(k—-DT'| 


t=(k-DT 


A derivada desempenha no caso de tempo contínuo 
um efeito de memória que a diferença desempenha 
no caso discreto. 
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Exemplo: Veículo espacial 


Lei de Newton: 


Força força = massa x aceleração 


(Sem atrito, simplificando) 


MO a ÇÃO, YO 
dt dt” dt dt” 
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u(t) = ma(t) =m 


u(t)=ma(t)=m y(t) 
Aproximação 


e SA da derivada 
dy(kT) oo Y(k + DT) — y(kT) pela diferença 


dt M ih de 1º ordem 
para a frente 


PyUT) + GAI) AT) 
dt? E 
HKD AGA yCK+DT) = (KT) 
T T 
T 


Fe 
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my(t) =u(t) Y() = = 


YV(k +97) -2y((k+ DT + y(kT) o u(KkT) 
fa o m 
V(k+29T)-2y((k+ DT + y(kT) = DAKT) 
m 
y(k+2)= uh) +2y(k+1)- (ke) 
ou 


VI = uk =) +2/(k=-D-y(k-2) 
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1.9 Conclusão 


> A noção de sistema permite exprimir pelo mesmo 
formalismo muito fenómenos de natureza e complexidade 
diversas. 


> Os sistemas têm em geral memória, que influencia a sua 
dinâmica (capacidade de evoluir). 


> A memória pode-se exprimir por diferenças temporais 
(tempo discreto) ou por derivadas e integrais (tempo 
contínuo). 
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Capítulo 2 


Modelização 
de sistemas por 
equações diferenciais 
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Vpa= CoaPpa 


TEC Ps Siste mas 
elétricos 
análogos 

Vaz Cura 
QR.= Pa— Ps Q R Je 


Fig. 5. Equivalent circuit for an aortic segment. 


https://snet.elth.pub.ro/snet2004/Cd/circ/circ P16.pdf 5/08/2022 
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ON, (1) & 


a b' 


(Witten (Ed), An electrical equivalent circuit model of glucose-insulin kinetics 
during intravenous glucose tolerance tests in dogs and in man, p. 1188)) 
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Modeling and Characterization of a Valved Glaucoma 
Drainage Device With Implications for Enhanced 
Therapeutic Efficacy 


Tingrui Pan*, Matthew S. Stay, Victor H. Barocas. J. David Brown. 
and Babak Ziaie 


IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, VOL. 52, NO. 5, MAY 2005 


Trabecular Coma 
Mesthwork (TM) 


Episcteeai mmina 4 Chamber (AC): 


Clliary 
Body (CB) 


Fig. 1. (aj A cross-sectional view of a human eveball with (b) a close-up view 
of the aqueous humor outflow pathway and circuit diagram representations of 
(c) normal aqueous humor outflow pathway. (dj GDD pathway, (e) an implanted 
GDD after fibrous capsule formation and, (f) a simple concept of bypassing the 
GDD. 


(DADC/DEI/FCTUC/LEB/MCPF/DPBF/2022/Cap. 2 


62 


R,, L, 
ES?) . o(9 
Ef) | 40 E 
RLC circuit 
A. B. 


Fig. 2. A - The electrical model of the respiratory system. Remarks: 
Eg(t) 1 Ext) — unknown excitation. generating flow and volume 
curves O(f) and V(t) respectively: B - Elementary section RLC used 
during the respiratory system synthesis 


Sistema elétrico análogo do sistema respiratório 
(B. Juroszek, J. Stanistawski, Synthesis of the Structure of the Respiratory System in Forced 
Expiration, MEASUREMENT SCIENCE REVIEW, Volume 3, Section 2, 2003). 
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Sistema mecânico 
análogo 


Modelo de uma perna (McMahon) 
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x(t) y(t) z(t) 


Intestino Sangue 


X(O) = 


Paes E 
na 


wv(t) z(t) 


Sistema fluídico análogo à absorção e excreção de um fármaco 
pelo corpo humano (de Bruce Biomedical Signal Processing and Signal Modelling, 
John Wiley and Sons, 2001). 
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2.1 Sistemas elétricos 


2.2 Sistemas mecânicos 


2.3 Sistemas térmicos 


2.4 Sistemas fluídicos 


2.5 Sistemas químicos 


2.6 Sistemas biológicos 


2.7 Sistemas fisiológicos 


2.8 Simulação de equações diferencias e de diferenças 
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2.1 Sistemas elétricos 


——— Componewks 


E e 
fifa q pe 
G— 17» —» << + 


. (a 
- tus : Vlt)= Rat) vtt) = e /ósar vrtt) = l id 
1(€) =C Ave) 
dt 


— parímetro: R,t,€ 
— armazenadara, Ads tverpi : condensadaras e bobinas, 
— apuanh: dad A evupio armazenada : 
condensador : he Cu? (enerpia bobnmil) 
bobina : % Lif Lene dinbhia) 
Se dissipodava, da 2 vai pra. : reg 
— grade: dada da upa dusspeado . . = Rise (miok 


Es 
e É ca metal par amslm de gilima: 
É bens ps A gpremer sersart pit ce correnlas (wos) 2 Unsses (malhas) 
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Exemplo 1 

SE qe 
ef () c po É 
A 


dúlt) a bo as o ii 
Cit) + LÁ, A fem at utel=o » RCy+LCy+y-Uso 


o | 
j+Ê y de ydo| 


Uma equação diferencial de 24 ordem 


ilt)= Cy 


Duas equações X=) > &=J="—YJ-—)J+>——U=-— 1X —-——& +— 
diferenciais de 1a L LC LC L 


ordem ps o 
EV e Ge Vves 
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Exemplo 2 


Tema dn wos 
0) - seleccons-e um nó cla 
referir (o mass, em prod) 
“x) — Oo densas mo eulro nas 
astnlhem-sa Lorvo Antéfiniba, 
CêL) - 4ShAVR-Se Jana Lgastço 


da KucheA (du corruda) 


Vo 


) (BD Lo ns du refentuur pasa tada ns 
v4) do ninisprrtos puros Cv) - eeprinmem-m as Corrudas 


V, 2 V; (1; =", 4 evlindo,) em fancos da, snipmis, 
2 dm frnsnelino da coxanito 
cus) Eguação no no (À) 


(e — Ap + 43 NÃ -me yo pl muro d Spa 
Equatp no no (8) da ns deve sa spas ao rúwero 
42 +4+44,=0 da smisprloo. 
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Exemplo 2 


Cv) UV = 3 a —+ C, Av; ouus dc / 1 fixo 
ds csoiecentas Ca ÁS 


E e A(E-w) 


RR E dig o 
R> 


df 
Escuvtudo an sv facão de UEM de Vo SW-V 
(agora = Vi-Vê,) 


Ve — Vez + Ve, -(W—Ve,) aih é lv, 
A, 2 
- (VU; — Ve 
“e (w é» ) C> d (VA Ve ) =. 
É) Ad 
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Exemplo 2 


E 4 RR | 
Escolta esa rage) 


VW o 1 vT 4 Pp ( N: 
REDE = Coe A 
Ca Cite “CR ss. «Crê 


N 


Duas equações diferenciais de 12 ordem 
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2.2.1 Sistemas mecânicos de translacção (movimento retilíneo) 


Ê a ste 


4 TI é , » E 7; 


atrito cutecca B mels ou 
cute «testuw dede ke. 
CNfm) 


Aplica-se uma Jota LH), jnodeg-* sam duditermudo gt). 


(em raposa, avtis de aplica LB), 0 dedicam ho 2 pulo e 
posar +" a de aajuímia v8o) 


- hun 
atak B male mete HM 
H(t) = 8 Sue) O=KIO fe)zM ad 


: pour nadam : B x; mM 
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2.2.1 Sistemas mecânicos de translação (movimento retilineo) 


- armazinadars du 4 do 


mala : =, K y* pf iantiad) 
mass : Pie tinta 
- dlissa padsr als 2 erga ã / ( eo ed 


atrito: w= By = B(g) (Mmjn om Wazts) 
(pinus) es/w.dk 


- prone fair fundo mute far a ansdcu 
x da Meto (7. da movimed ) 


soma de toda, as UAM Acluusemo à 
? Miva 


(1) (eg) Cm /p*) 
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Exemplo 1 


Quar ds friso aflicade, o mana [2 
eu, pousa de eetrrior 
o — alado opsê-s às ymavi- 
estado, due Sa, aflia AA VAR 
forca da sinal fndirdro asd U 
h = 8; 


ecko Mo AA. 
dec da Malan 


u-Buy =H4y 


.. s 
My +68 Y=4h 
Soma das forças aplicadas a M 
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Exemplo 2 


Quo as foras oflcadm a hº 
e mu, Apnicada da gelario 


o ul, += By 
e O sê exusa à muda sob 
Mah, fis 4 


ln da Newton 
U- Ig - 6y = ny 


M y+ by “ay = dA. Soma das forças aplicadas a M 


(DADC/DEI/FCTUC/LEB/MCPF/DPBF/2022/Cap. 2 


Exemplo 3 


* Suldensar da Mm cudbineval: 
à AM. M2 mass ontomevel 


ut forca velicol aflora as 
vultnd (po acos ola prsmioleda 
Anda da 


My +By+ky = AA 


(semelhade as teamo 2.) o Br ka, dl 
PRaita/m 


(2: drf. Crua, 2ordam, Cof. cont.) 
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i bj 
Cf tt 


fest) f=lh-%) “ul | 


12 fposicas ; y 2 velidada yr; y £ E acotorack q ; 


a=u 


t= ly) fé) f=hd 
! = K(Mj-4e) 


Re a: Ed 
8 x v M 
“ Va É 


Vi-Va é Ps o 
a Bl asp dE=pÊ fendronç 
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Analogias elétrico-mecânico 


Quadro da Amailopia escimito eo ele? 
Traholees 
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2.2.2 Sistemas mecânicos de rotação (movimento angular) 


— compontrias = Povo 
da 
B Tr de Ea 
atirck ou PR E, Ba nerio 
frícess fotebive 8 redes (N-m) 
(N-mfreafa ) 


Apluca - at uam bina T(t), produg-a um rota (ole sta - 
perda on pls ) O" em atlacg " pose Arjualas da 
£,2 


friceê Rê 


Corças kk 
Ttt)= 6 dot) . 
dk 


aneraa J 
Tt)=k0/t) Tt)=7 dot 
dt 


= B w(t) = Jx(t) 
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2.2.2 Sistemas mecânicos de rotação (movimento angular) 


- parônide U B, k / J 

-armezenadars de tuas: 
torças : = k e” (pinus) 
pvsDaR: O =1h Tu? = T(0)* (ont) 

— dussupado de nar A o sceç 
putcego pfaltuna nd =8 19'= B(6* (Mixia) 


E prewepar Arado veado de amsd.na dE” próiicd 
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K 2 cofiwul dito de mio 
BS juca evita o dia 
La Sajutice exterior 
T 2 binam'o estaria afhicade 
TS moment de cneru? em 
rela Oo Xe. 


Gunass 26X O) binanio aflitado do diste? 
o binosio exterie T(E) 
e benaluo resultam da torcãr da um T = 6tt) 
o bina davide as olaik B, To = Bélt)-Buwtt) 
dor da Nimilm aflicada as duste 


[dA Las J6 


O 


pisada — 76+R6+k0 =T 
9 + É O+ = g =L ; nda du. 2“*ndum, Lya anverianh 
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2.3 Sistemas térmicos 


Os sistemas térmicos transmitem e armazenam energia térmica 
(quantidade de calor) 


Exemplo: termómetro de mercúrio 


tmb: Tm dps à Wimfrtua Ta 
Comenta esto ds ojusa Co rt 


quantada da pala nominikiola aro neaatuinio sas 


tada snÃamiT : ' -— 
440) = Te Tonteo PER 
EG uaRrT N prada do tirinal, 
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Exemplo: termómetro de mercúrio 


Cutas mewk als E mjurduss do nustuno 
Timo) E Em puitus tuts do merulrio 
Tm(tl — Tmlo Exuurmenk AQ Enpitua 
C E csponwdd tirniw de munio (tslal) 
C-Tmlt) E quedas d tolm contula no mmairoo 
C (Tm (6) -Tmls)) € Quant: drolt du fede ebsom- 
du pe mutus toda O ol É 
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Exemplo: termómetro de mercúrio 


balanto de Lu pi 
qe). ot = a[cTm lt] = c.smét) 
AnonÃiAdada Ras ia ts 


To Tm -C O Tm à =- RC To tm =Ta 
e ot 4t —>o 


Uma equação diferencial de 1º ordem 
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Exemplo: cilindro doméstico 
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Exemplo: cilindro doméstico 

(O 9, E quaniidad du coler forrucido pala renatiuw uai 
O WE quantidade ole ester absmnda pula apar no colidir 
(D 208 Quantidade de ralor qusa: de utiudar (va aíjua) 
(O q, 4 quantidade de cale qu sudás pa lindo (na ala) 
O dy E Qunelidad de ada. apos forma (pula pride) fu fon. 


PRINCÍPIO DA CONCERVACAS DA ENERGIA 


fe pio * Mm Gt Gota 
calor qua auitve dias cola Gal Sãn 
o fred 
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Exemplo: cilindro doméstico 


F 4, : frneudo pula resina aladanis 
em cada «valawia, 
dp = Ra 2 4u 


- Gp: asno pela pus vor rextindom 

durant um antinvaho dá nho AE 
4, = M.c.oT 

mam ds = 
apra tSpuiitic dep (Qua. da altus 
EA - 
ty = Ee: Ta Neo E YU - 

dim 6t-»o 


H AT 
“de 
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Exemplo: cilindro doméstico 


- 4: metido na apar apar euntis com Comet 
, ep A LÊ 
Le = Vel 


-9: AGA LOM. & apa, de Ma tundro 
Jo -=V.eT (adwiliudo pu 0 Condo 
dt nad e'sprslas de eubiado) 
—-9* prde-ne fps o omni alias da fas 
R$ rsdinve term da hauda ThTa 
TR %a 
Gas R 
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Exemplo: cilindro doméstico 


subtulwinde na E do bolanto de trupe 


ur VeT = me Lev. eT + EE 
Pald quo Ape T VT 


Uma equação diferencial de 1º ordem 
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Analogia elétrico-térmico 


Resdinus Lrmia 
4 = LA(T-Te) 
q = : (=) 


A 
Tt ' R$ 4- ai al Te 
-. * 


— 
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Analogia elétrico-térmico 


Termito de= Elcetrio 
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2.4 Sistemas fluídicos 


Nos sistemas fluídicos circulam fluidos, ou seja líquidos os gases. 
Neste curso interessam-nos os líquidos. 


, — mt) 
A, Aa = seccoe do Veor 


ê Pluvolicas da, vellvulo, 


q. 


Pal prassegum 
Ep %e 492,95 & como, 
k ) pa 
is Ss Eos Aula 2 mis da bd 
' ê 


Ay BE) = Gult- go6e) 
ArÊaMe) = 946) -92(E) 


Ver os apontamentos, págs 46-48. 
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Sistemas análogos: sistemas de natureza energética 

diferente (neste caso elétrica e térmica) mas com a mesma 
estrutura (componentes que têm funções análogas e sinais 
análogos). As suas equações diferenciais têm a mesma forma. 


Dois sinais são análogos se, embora de natureza diferente 
(corrente e calor, por exemplo), têm a mesma evolução 
temporal quando os sistemas têm a entrada com a mesma 
forma no respetivo tipo de energia. 


Dois componentes são análogos se exercem sobre os 
respetivos sinais de entrada o mesmo efeito. 
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TABELA DE ANALOGIAS 


2.5 Sistemas químicos 


= Envolve os princípios básicos de termodinâmica, cinética, fenómenos 
de transporte, etc. 

«= Um sistema químico recebe materiais e energia do exterior, fornece 
materiais e energia ao exterior e dentro dele acontecem reações 
químicas. 

= Aplica-se o princípio geral de conservação de massa e de energia 
(nada se perde, tudo se transforma). 

« Genericamente, o princípio de conservação do componente S do 
sistema é dada pela equação de balanço. 


a de S | ps de S para | py de S dio] ip de S | pao de S | 


dentro do sistema dentro do sistema fora do sistema gerada no sistema consumida no sistema 


ES ss aÃ ASAE MMaSEJSã Ss 
intervalo intervalo intervalo intervalo intervalo 
de tempo de tempo de tempo de tempo de tempo 


S pode ser a massa total, a massa de um componente individual, a energia 
total. 
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Assim, para um sistema químico genérico, poderemos escrever 


Balanço mássico total: 


MD. ano o 


i.entradas j:saídas 


Balanço mássico para o componente A: 


d(n,) d(C,V) 
e — o = >. Cas > Ca FP; trV 


i 
i:entradas j:saídas 


Balanço energético total: 
dE d — +K+ P) o 


Ei > oEh- > pFh +Q+W, 


iientradas j:saídas 
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p & densidade do material do sistema 

p; & densidade do material do fluxo de entrada i 

P; à densidade do material do fluxo de saída j 

V à volume total do sistema 

F, £ caudal volumétrico da entrada i 

F, à caudal volumétrico da saída j 

n, É número de moles do componente A no sistema 

c, & concentração molar (moles/volume) de A no sistema 


A m . 
c, = concentração molar de A na entrada i 


A 
CA, à concentração molar de A na saída j 

r É taxa de reação por unidade de volume para o componente A no sistema 

h  entalpia específica do material na entrada i 

h, à entalpia específica do material na saída j 

U,K,P& energia interna(U), cinética(K) e potencial (P) 

Q é quantidade de energia trocada entre o sistema e o seu ambiente por unidade de tempo 


W, & trabalho (shaft work) trocado entre o sistema e o seu ambiente por unidade de tempo 


Aplicando a um sistema qualquer obtêm-se equações diferenciais de 1a 
ordem, que constituem o seu modelo matemático. 
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Exemplo: aquecimento de um tanque bem agitado (stirred-tank) 


Seja um tanque que está a ser alimentado por um caudal de líquido Fi (m3/min) 
à temperatura Ti (ºC). No tanque é aquecida por uma serpentina a vapor onde 
entra um caudal de vapor de Fv kg/min. Do tanque sai um caudal F de líquido à 
temperatura T. O tanque está bem agitado e por isso a temperatura da saída é 
igual à temperatura dentro do tanque. O vapor de aquecimento, ao perder 
energia na serpentina condensa, saindo o condensado para o exterior. 


hei dh F Ff 
dt A A 
dT l Q 
o a ( Lo ) + — —— 
dt Ah Cc 


Note-se que este sistema 
químico é a conjugação de 
um fluídico com um térmico 


"ot (de Stephanopoulus) 
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2.6 Sistemas biológicos: Modelo Lotka-Volterra do par 
predador-presa 


Os sistemas biológicos descrevem o funcionamento dos seres vivos. 


Este modelo ( https:/Awww.math.duke.edu/education/ccp/materials/engin/predprey/pred2.html 5/08/2022) 
descreve a interação biológica entre duas espécies, na relação particular em que 
uma, predadora, come a outra, a presa. Há muitos exemplos na natureza: leão- 
gazela, pássaro-inseto, pandas-eucaliptos, etc.. 


Hipóteses simplificadoras: 

> a espécie predadora é totalmente dependente da presa (é o seu único alimento). 
> a presa tem uma quantidade ilimitada de alimentos 

> a presa só tem um predador, o que está a ser considerado e não tem outras 


limitações. 


Na natureza as coisas são mais complicadas, mas assim poderemos desenvolver 
um modelo aceitavelmente complexo. 
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Se existirem zero predadores, a presa não tem qualquer limite ao 
seu crescimento e por isso cresce exponencialmente. 


Sendo x(t) a população atual no instante t, a sua taxa de 
crescimento será, sendo a uma constante: 


axo) 


t = ax(t) 


Mas havendo predadores, o crescimento será menor: os 
predadores introduzem um factor negativo nesta equação. 
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Se y(t) for a população de predadores no instante t, cada um 
deles tem certa probabilidade de encontrar uma presa. Essa 
probabilidade é tanto maior quanto mais predadores existirem 
e tanto maior quanto mais presas existirem. Por isso se 
pressupõe que 


> a taxa de encontros entre um predador e uma presa é 
conjuntamente proporcional ao tamanho das duas populações 


> dos encontros realizados, uma parte b fixa resulta na morte 
da presa 


AXO xt) bx(t)y(t) 
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Quanto à população de predadores, y, considerando-a 
isoladamente, não tem alimentos (não há presas) e por isso 
morre a uma velocidade proporcional ao seu tamanho, ou seja, 


neste caso, 
dy(t) 
= = cy(t 
+ v(t) 


por não haver energia para novos nascimentos. 
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Ora felizmente para os predadores existem presas, e por 
isso eles vão ai buscar energia para se reproduzirem. Por 
cada encontro predador - presa há uma parte p que 
resulta em alimentos para o predador e por isso a sua 
população variará de acordo com 


E =-cy(t) + px(O (1) 
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Juntando as equações obtidas obtém-se o modelo de Lotka- 
Volterra 


a = ax(t) —bx() y() 
mo = -cy(b) + px(Dy(t) 


composto por duas equações diferenciais de 13 ordem 
acopladas, que não podem ser resolvidas separadamente e 
para as quais não é possivel encontrar uma solução analítica 
(mas apenas computacional). 
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Pode-se mostrar que existe um par (x., Yy.) para o qual 
dx/dt=dy/dt=0, isto é, as populações estão em equilíbrio e 
não variam ao longo do tempo. 


Para isso basta fazer 


O = ax(t) — bx(tD) y(t) 
O = —cy(t) + px(t)y(t) 


I la-by(DIx(t) = 04 x(t) 0 | [a-by(t)]=0 VCD) = 
=> —s =" 
[-c+ px(OD=0n7(0%0  Il-c+px(t)]=0  ug- 


sn cla 


um tal estado chama-se regime permanente ou estado 


estacionário. O seu valor depende das constantes a, b, c e p, 


os parâmetros do modelo. 
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Existem outros modelos de duas ou mais populações. Em 

https://www.montana.edu/screel/teaching/bioe-370/documents/L-V competition 
populus.pdf 5/08/2022 

encontra-se um modelo de duas populações competindo no 


mesmo ambiente, (a que os autores também chamam de 
Lotka-Volerra), mas neste caso não se trata de relação 


predador-presa. E em 
https://pt.scribd.com/document/498569948/25-MultiSpecies 5/08/2022 
um modelo de várias espécies de peixes competindo pelos 


mesmos recursos. 
Existem outros modelos que exprimem a: 


> interação agressiva (combatem-se ativamente, embora não 
se comam), 

> cooperação ou simbiose (ajudam-se mutuamente), 

> relação predador-presa, 

> interação forte-fraco (uma espécie está mais capacitada 
para SORLev INE), 


Ver DP. EX. i/Biological interaction 5/08/2022 


2.7 Sistemas fisiológicos 


Os sistemas fisiológicos compõem-se em grande medida de 
sistemas dos tipos que vimos (mecânicos, fluídicos, químicos, 
etc.) que ao incluírem órgãos vivos adquirem, por sinergia, 
novas propriedades. 


Usam-se frequentemente as analogias com os circuitos 
elétricos para estudar as suas propriedades. 
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Exemplo: um vaso sanguíneo 


V É volume do vaso | = 
Propriedades essenciais: 
Q 2 caudal de entrada 


PAN my 1 
P = Pressão à entrada do vaso Resistência : oposição à 
A ; circulação 
Q, = caudal de saída 
É piessisd saida Complacência: aumento do 
2 volume com aumento 
e é A RN : 
P (Pressão exterior) = 0 (referência) do caudal 
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Resistência do vaso 
p = () ( de Hoppensteadt,9) 


Q, (Q, 


o P, V P, 


Para um vaso sem complacência, rígido, de volume constante e bem 
conhecido, teremos, como nos sistemas fluídicos, 


2 H-b 


“= R 


A R chama-se resistência do vaso e teremos aqui um vaso de 
resistência. 


(DADC/DEI/FCTUC/LEB/MCPF/DPBF/2022/Cap. 2 


Complacência do vaso 


Num vaso elástico sem qualquer resistência fluídica, a pressão na 
entrada é igual à pressão na saída, para qualquer valor de Q, 


P=P,=P 


A relação entre a pressão P e o volume V do vaso pode ser aproximada 
pela relação 


V-=CP 5 Cc! 
P 


em que a constante C se chama a complacência do vaso. 
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>» 


Q 
SRA, 
eita 


Um vaso rígido, puramente resistivo 


= pm - -—— (o — 
-— 


» 
£ 


Um vaso puramente elástico, sem resistência 
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complauúununa € Capruidor CO 


V=CP Q=CUV 
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VIZ 


Resistência e Complacência simultaneamente 
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Q; Q, 


= dy, NV, 
P, variável Ca O 


Vaso Circuito eléctrico 


Complacência, C C, Condensador C 
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2.8. Simulação das equações diferencias 


Três blocm deutiinica betviess : 


Ur “ Ve 
[A 
Va 


"; VAU; vo =X. di ta. 
Semader Am blkivador grito "ua 


o je JA > s 
Gain Integrator 


Corresponding Simulink blocks 
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Exemplo: Equação diferencial de 13 ordem 


jrey=u O 


j=-2y A 


(ea. dg. da 1º oroum ) 


eg di. du (É avelma — fw nessas um arliarado 
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Exemplo: 2 equações diferenciais de 13 ordem 


1+Sp=u 4=-SY +4 


E-3Y eso 
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FIZ 


Exemplo: Equação diferencial de 23 ordem 


J+95f-29=4 SI, 


Y=-g5y+2g+4 y=/qdt= [q dba 
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Exemplo: 2 equações diferenciais de 22 ordem 


ray Bj ros n+0N% = 4, D= 29, 839084028 +4, 
da +91 =) “Sp +ISh =2u, U=-915,+9, +58, -48 hp +20 
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Simulação de equações diferenciais não lineares 


Modelo de Lotka-Volterra Diagrama no SIMULI NK 


O ax(t)-bx(t) (O) 
dt 
a = -—cy(O + px(O) y(t) Osciloscópio 
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Simulação de equações de diferenças não-lineares 


Modelo do crescimento de uma população isolada 


Xe — Ax Ux) = Ax, o 


Diagrama no SI MULI NK 


Visualização 
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Conclusão 


> Na generalidade dos sistemas podem-se aplicar os 
princípios básicos (do respetivo domínio) para se obter 
um conjunto de equações diferenciais (nos sistemas 
continuos) ou de diferenças (nos sistemas discretos), 
que são uma representação matemática do sistema e 
constituem por isso um modelo matemático. 


> Esse modelo servirá para o estudo do comportamento 
do sistema (a sua evolução temporal, por exemplo) e 
contém caracteristicas matemáticas que exprimem 
propriedades internas essenciais do sistema. 


> Obtido o modelo, importa saber como e para que se 
vai usar, e é esse o tema dos próximos capítulos. 
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CAPÍTULO 3 


DINÂMICA DOS SISTEMAS 
BIOLÓGICOS E FISIOLÓGICOS 


Do capítulo 2... 


Lo » | rm) fra =!) a 
a” 2 UR + AJ HAI = DT) + bay +. tby4 +b, 4 


rr yo | 
Ala com do same pulicda u qual, Gun o Y oblidl 
Lo resobucal da qua dofranial 


.. COMO ??? 
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TRANSfSR ADA DE cneence do 
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3A.4. Transformada de Laplace de funções elementares 


Impulso de Dirac 1 
(1 [,0,t=0 I stjat=" 
O )=s e lo É = 
( 0,t=0 sá 


h 


to) = É [66] a [ Cote) Et = 2º] ma 


[7] = 
(tod. o E) sd duto dn lipha ol cube psem) 
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i) Degrau unitário 
l 
(o) / 1. t>0 
[j= 
As lo, t<0 


ato) = £ [att] = [Co) Cat a [ tab 


E *º| E -5[º = | - 
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li) Exponencial 
o cal 
Exponencial L | 


ptiy= LD] = [Perde = [Ea 


= 4 ea im - sá [e-1 


H" 


pa A+ a 


o 


/ 


OO 


A+a 
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iii) Rampa unitária 


anda mto Si)= E 


Ft) 


(era (o 


4 
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3A.4. Propriedades da transformada de Laplace 


1) Linearidade 


filo): É (464) + lo) = 2 (e,(41] 
? [e4,()+8 bt4)] - “hlo)+ Bolo) up das 


ii) derivação no domínio temporal 


t[ <ue)=2[ito)=2(e8)-ge) = 40). ph) 
dio fu] =4Fln)-» VC Elo)-n' Ho) -n Lo)... Lo) 


male  flt|=o L£<o 


din we devultino Wanpasl es nudh ta 5 
Em ts A emp ar eai 
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iii) integração no domínio temporal 


Lt pato ne domino Abdinal — dota pra 
br domino Lehltio 


A [ / Ztejár] = L Ela) 
pe 


(DADC/DEI/FCTUC/LEB/MCPF/DSBF/2022/Cap. 3 132 


iv) deslocamento no domínio complexo 


Peslecameks wr da yum Corafrliso, ac prada hdoes 


f [2*pte))- F(n+a) 


viii) deslocamento no domínio temporal 


eine no dana inio dive funal — mualALcaeg 
fm & A UA demands Compltxa 


É [flé-a) = 2" F(4) 
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iv) Teorema do valor final 
Este teorema relaciona no comportamento de f(t) no infinito com o comportamento 
de F(s) na origem 


Eae fÉ) = bruna A Ela) 


Do (24 Fila) nês tm pelo mr SPD) 


v) Teorema do valor inicial 
Relaciona o comportamento de f(t) na origem com o comportamento de F(s) no 
infinito. 


Lima He) EE po we Afla) 
Naa (onto) 4 Gloom foi? 
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3.5. Inversão da Transformada de Laplace: 
decomposição em frações parciais 


Uma transformada de Laplace é em geral uma função racional, dada por 
D(o) ua): pdindrdr de o de pn 
D(A) : fretimas nem dA da “4 


Diz-se que ela é: 


Ri) testei 7 gras Áflo) € you pd e 

as Aa mbinÃT 2: v ta) € D/a 
Fo) bibréproa : o dio a dá asa 
Ha) atadaRlhria, : Anata Ma) > Graus D/o ) 
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Definição de pólos e zeros fe F(s) 


(A) = AACIA 
D (0) 


N(s) e D(s) são polinômios co-primos (não têm raizes comuns) com coeficientes reais. 


dell. ponto da epuaeso D(a) = 


* VOL da epuata Nfn)=e 
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Como sabemos. uma equação de coeficientes reais pode ter raizes 
* reais simples 
e reais múltiplas 
* complexas aos pares conjugados (conjugadas) 


e por 1sso os pólos e zeros podem ter qualquer dessas formas. 


Quando F(s) é própria. sendo N(s) e D(s) co-primos. o 


au jorlem] ou FG): mimero dt bolos 
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Exemplo: 


ê 
El) = O antenas | bihrobres 
At + TA + IE NEON +$ 


4) furtarigao de dimencinaaa 
den =[14,%,18, Zo, 8], ? = roots (den) Matlab 
cá P 2 [4,2,2,2) 


nt, Va HIS 4204 48 = (p en)(n42) 


| Eu pl rd praalb:) eridada 3 
“fail? 
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UA) dacomposicas sm fracçeio qo 
Fada A — pda dy So q ha. 
At pe =, 


bedo sÃão Lima perdi ma tabela 
44) snversas 
fU= A EUA, E a ha DAtE, Ag EE 


Cocina 02d tada hs, A, 2, Ps, Pas exllonk ela Nécôlms 


(DADC/DEI/FCTUC/LEB/MCPF/DSBF/2022/Cap. 3 139 


3.6. Aplicação da Transformada de Laplace à 
resolução de equações diferenciais 


Seja por exemplo a equação diferencial de “ordem 


ay +by+cy= A 
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Aplicação de L a ambos os lados da equação: 
E [ay +by +Cy]= 2[u] 
cx) fala Lovatdeda da É 
A ÉTy]+b E[9] +e Ely) = 2 [u]) 
uy ju drama da dica, 
ET): A Liv) to) - 906) = 0º Yo) aylo)-jfo) 


Ele 2bfy]-gto) =AYHo-ybo) 
ZIy)= Y4) 
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ay +bo+Cy = 


U) 


nº a Yo) - na yo) -agto) + ab Flo) - bah) 40 Ylo) =U(0) 
(a at sb» +c)tilo) = (ap +b )ylo) + a gfo) + U/») 
Tn) = (0 0+)90 +ayto) |U) 


z | 
antebnra. ap bapo 
sob olevida às Cmdia panda cenda 
Luas , Var à dando, Yap 
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Tabela das Propriedades da Transformada de Laplace 


Property Time function Laplace transform 
Linearity o fi(t) + Oo fa(t) o Fi(s) + oo F>(s) 
s-shifting eEf(t) Fís +a) 
Time shifting f4-NDgut-T)forTz0O esTF(s) 
s-differentiation (e) SE (s) 
Time differentiation dep (1) sF(s)- f(0-—) 
4 
dº > = p 
ae 7 (0) s F(s)—-sf(O )— f(0) 
Time integration In f(mdr F(s) 
0 s 
As 
Time scaling flat), fora >0 PR 
t 
Convolution Í h(t— Tu(T)dT H(s)U(s) 
q 
Final value f(co)=limsF(s) (if sF(s) has no open RHP 
pás poles) 
Initial value f(0+) = lim sF(s) (if F(s) is strictly proper) 
S—309 
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f(D,120 
o(1) 
lorg(t) 


t 
t” (n: positive integer) 
e“ (a: real or complex) 


te 


id Mud 


F(s) 


que 
pr 
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Tabela da Transformada 


Sem oscilações 


sin Dot 

cos Mot 

ft sin Qt 

ft cos Dot 
e* sin wgt 


e“ cos mt 


y 


om oscilações 


Do 


s? + 6 


Ss 


s2 + o 


2Wos 


(sé + og) 


2 2 
s* — Wo 


(s2 + wo) 
Wo 


(s + a) + mo 


s+a 


(s + a) + w5 
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Exemplo (Bruce): 


Imagine-se um músculo de rã fixo por uma extremidade e no qual se suspende 
uma corpo de massa M. O músculo, devido à sua elasticidade, alonga-se. Pode-se 
desenhar um equivalente mecânico de translação com uma mola K (elasticidade), 
atrito viscoso B, massa M. 


A força aplicada à massa é a da gravidade e igual ao seu peso: F=Mg. 
O atrito e a mola reagem com as forças 


dy 
atrito: F.=B- =B 
B di y 


mola: EF. =Ky 


(b) 


A força resultante é a soma das três e 


produz o deslocamento y. Pela Lei de muscle ; 
Newton : B 
Mg-Ky-By=My >My+By+Ky=Mg 
co e weight 
mm (mass = M) 
(Fig. 2.10 Bruce) Mg 


Uma equação diferencial de segunda ordem, linear. 
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n» " =M 
My + By +Ky=4 lido 
Se, para fins ilustrativos, supusermos M=2, B=3, K=1, vem 


29434 +Y=U qué a(hl= + 


Nota: se se considerar g=9,8ms?2 , então u=2*9,8=19,6, ou seja, uma degrau de amplitude 19,6. No entanto um degrau 
unitário é representativo para a resposta do sistema; é como se o caso se desse num planeta de gravidade g=0,5ms?. 
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Aplicando a Transformada de Laplace 


2 (2º Yto)- 1 yo) -glo)]+3 LoYto)-yt0)]+ Yin) = Vo) 
(20º+30 +90) Yo) = 2a) + 29/60) 4396) + Ulo) 
Yo) = (easy) tado, A ut) 


ZA? +36 +4 2843044 
siri E 
dada: Flol=-4 * Ru fas 
Glo)=1 
(2043) +2 
O ——— O ————— 
20º+30-+4 “(24 “er + TED 
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ZS 


Z 


7 
Ea 
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2 
Polo ESA, 


A (041) [0+485) 
ah, A 
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3.7 Polinómio característico, modos do sistema 
e estabilidade 


Se entrada nula, Y = Y,, resposta não-forçada ou resposta natural 


vio) = (2242790) tz ito) Yo. TOR Ãa. 


ZA" +30 +14 +85 A+i 


Modos (característicos) sto + 


ost -É Ja (tl) = As L i + Az L 
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A resposta natural é uma combinação linear dos modos 


Tec) = 2fato=g69] 0 [-gto) ago] E 


(DADC/DEI/FCTUC/LEB/MCPF/DSBF/2022/Cap. 3 152 


Estabilidade (resposta natural) 


Se um dos modos tender para 


) A, Aa. infinito, a resposta natural 
Tas (4 = gi tende para infinito, o sistema 
n+0,5 A+1 diz-se instável 

Condição necessária e 
suficiente de estabilidade: 
-0st A - Seas raízes do polinómio 
y (&) 2 AL " + AL forem todas menores que 
A zero 


* Seas raízes do polinómio 
forem menores que zero 
excepto uma (e só uma) 
que pode ser zero para O 
caso da resposta natural, 
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Resposta forçada (condições iniciais nulas) 


Yo) = —— Ulm) G (0) Ula) 


ZA +30+1 ia 
E 4 
fz4 FUNC DE TRIM FERÊNCIA 


A (2) = a 


(a+40 8) (n+1) 
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G (o): 


(a40, 2) (n+1) 
Ulo)= Elgin] = (oupura amltaiio 


Vala 6 (o). Vl») = 


Ta a 


Í 4 A3 
zo +4+—— 4 —— =.24 Ê | | , 
P +49; 9+1 9 PAGE pl : 
id E + 2 ai” /Êm0 
o dani om dei do de sm 
lado, NV; (n ) 
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Estabilidade (resposta forçada) 
gltl= tr a sa sHÊ tro 
E E aa 
arado U(n) 
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Polo real simples 


Vais dos cos qul=LE 


Thls atos 4 gt = 2.58 


Llo) = 4 me 4(t) = Na 


Yo) x à. di qt) e | 
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Polo real simples 


E 
Lim yéla 6 


É -s00 


| F, 4 q 
Í V 
& + a 
tim q (| «00 
té -s00 


E E O 
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Polo real duplo 


1: + E Jtt)= to 


te 
(A 41) | gi Lat 
Ta) = er ad as Ná + 2 3 q 
+ 
Ta) = re 7 —S 96) = Ee* mn jo pit 
aÃ A 
nd e E z 
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Polo real duplo 
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Polos na origem 
Pólo simples em G(s), entrada em degrau 


| o. 
Y(s)=>>y(D=t>limy(t)=0 4 
a Instáveltt 
Pólo duplo em G(s), entrada em degrau | | 


Y(s) = - > y(t) = -É => lim y(t) = 00 


t—s00 


.. Origem é um ponto de instabilidade !!! 
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Polo imaginário simples 


Tide Lo. —s 6) nda f 
A +41 


XY) = Lo + yt)=fmat 
A 4 4 
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Polo imaginário múltiplo 


Ylo)= L2. — yt)= €4ent 
(15 cresunta com É 
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Um par simples de polos complexos conjugados 


Yo) = Ra + JU) = Ltsent 


p=-4 a 
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Um par simples de polos complexos conjugados 
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4 


E fa E 


" 


Resumo 
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Os zeros da função de transferência 


LA à UA E” 
e a ue) EE (2 + ) 


Pl AR o yti= [ecai ano] 


Ata)ln+4b) 
o zero altua cx totftends 
da, 2 honenuar, 
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Conclusão 
Jtt) é uma combrnacas bineea e pepoventizis 


os polos dofinum à ma dao teorema 
é (acças dalina ) 


/ 
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3.9. RESPOSTA TEMPORAL DE SISTEMAS LINEARES 


Sistemas de 1º ordem 


Sistemas de 2º ordem 


Sistemas de ordem superior 


(Nota: pressuporemos que são todos estáveis) 
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Sistemas de 1º ordem com entrada em degrau 


Glo)a —K cm 
ALA 


A+a 
U(s) == (aso) 


=== Fa? É (1-2 = 


A 
Plana s valer final (em Taco) 
qu - 
a ni A | Saída y(t) 
5 a 
canadindicas 


Ra du têmpo te 
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Sistemas de 2º ordem com entrada em degrau 


Saída y(t) 2? 


2% >1 >» & polos rear Simjls, 
4 4 =4 = 4 bol, pral dino 
ocê <1 id Aonde = 


Aya a -$y ty 
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Saída y(t) 
4) 421, 2 bom reais distals 


Açyê — 2 A pe wu V$2-1 


Sé | 
AFE3A+2 Plava 1 
o E 
ak Ss 
P+D(0+2) o; 
se 
s para diaÃd a 
J Gare K/2 
e qelre: dasta as subiola, R firma da pum 


(depusndim du 8 e tn). 
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vu) G=2, um pelo acl dufho 
Exa 
Gn) = K KR 


At+20+1 (a x 4)€ 


ye 
Plana k 
a 
” a 
a o 
(0 
(08) 


caraÃipiÃ an 
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U(s) -— 


LA4) o<4 ct, 2 polo conflixao conjugado 


sk 
Gl asa een pi 
AP +27 A+ tuR as 


E EE 
pragas 


— t 
PO O pe 
POE 


= castf 
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U(s) -— 


Pra == Um tj omnVi- 5º 


| fe o. 


qto) = (o) 
/ Y(c0) = o? 
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É 
G(n)= a 
AÉ+42S UA + i,d 


Saída y(t) 


Oder: resposta sub amertêusa 
É 1 » nechordo om amntêcmenk Cr chiur 
291 : rshedz sabre -amnTacts 
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2 
G(») "S.A U(s)= 1 
S 


+ + k 


T= un) constant du tomp PO =100«p(Tl. ) sobre eliane em / 
To = 3T (st) Aêmpo cu estabelec 


no ” - 
Te seg. emp da putida kr od 
To: = — tim bo ala pros 
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Sistemas de ordem superior com entrada em degrau 


pigs 
4 Polo! ADC oa di! ; 
PO — 


Saída y(t) 


T(o+3) 


ai (046) Yo+b (espe Josh) 


(Nota: sistemas estáveis) 
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Álgebra de blocos (132) 
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Capítulo 4 


ESPAÇO DE ESTADOS 


Exemplo 1 


No modelo do músculo da rã suspenso (Cap. 3) obteve-se a 


equação diferencial de 2º ordem 


My+By+ky=u 
Pode-se reduzir a duas de 1º ordem, por 


substituição de variáveis: 


= 


MA = do 
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weight 
(mass = M) 


(Bruce) 


N o o 
KM =>YDX=Y=»* 


p º o o co B e K 
M=yPM=0)=y=="",)-0) +04 


M M M 
á B K ] 
OU sgja X=-— X-— xXx +—U 
M M M 
A = & 
' B 
dd x +—u 
0) ] O 
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O que queremos observar ? 
A posição y ? 
Neste caso a saída é y e portanto 


Y=*, =1x,+0x,+0u=[1 0)x+[0Ju 


y = Cx+ Du 
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Exemplo 2 


No modelo de inspiração-expiração pulmonar 


V(t)+aV (t) = bP(t) 


O que queremos observar ? 


V? Nesse caso x é a saída 


xº2V uêP(t) y=* 


x=-—qx+bu y=x+0u 
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Exemplo 3 Metabolismo da glucose 


ne —-mg-m,h+J(t) 
a -mh+m,g +K(t) 


g: desvio do nível da glucose do seu valor recomendado 
h: desvio do nível da insulina do seu valor recomendado 
J: taxa experimental de infusão de glucose 

K: taxa experimental de infusão de insulina 


m,,m,,m,,m,, constantes características de cada indivíduo. 
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Estados: xãg  xêh 


Entradas: u é J(t) u, é K(t) 


Saídas:  Yy=X o )=»% 
dx, 
=MMX-—MAX+uU, Y, = 4 
dt 
dx, o o 
E =+m,X—m,x, +U, XY» =X; 
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Exemplo 4 No modelo de Lotka-Volterra 


a APP E O 
CO ey()+ Ovo 


x, É x(presas) x, 2 y(predadores) 


dx (t) 


= ax(t)— bx (Dx, (t) 
dx(t) = -—cax(t) + px (tDx(t) 
dt 
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X 


l 


Z 


= 0X, — DX,X, de (XX) 


= CX, + PXX, E? = 5(x,,%) 


Sistema não linear autónomo (sem entrada exógena) 
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Caso geral 
e fi QD, cs (O, U(D, o» Un(D, O) 


YO = 9 CO, o AD, UMD, +. Um O, O 


n variáveis de estado 


RR 


f(x,u,t), x(to) = x0 
g(x,u,t) 


Caso de sistema invariante: o t não aparece explicitamente 


a o 
| 


bs 
| 
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Equação de estado x = f(x,u,t), x(t;)) = x0 


Caso linear invariante X = Ax+ Bu, x(to) = x0 


Equação de saída 


y = g(x,u) 


Caso linear 


y = Cx+ Du 
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= Ax+Bu, x(t;) = x0 
= (Cx+ Du 


Condições 
iniciais 
(memória) 
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x2 


Plano de fase e curvas de fase 


x1'=-2x1+3x2+u 
x2'=-0.5x1+2x2 


Plano de fase e curvas de fase 


x1'=-2x1+3x2+u u=1 
x2'=-0.5x1+2x2 


traçado com o pplane.jar disponível em https:/Awww.cs.unm.edu/-joel/dfield/ 
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Resolução da equação de estado 


- pela transformada de Laplace 


— £-Ax+Bu Z [). A 2fi+8P[u] 
AXbU)-Xo) = AX(n)+B ULlA) 


— Ve Cx+Du LTy]= ctl]+DÊ[u] 
Ylo)= C xls) +D Uln) 
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A Klan) - X lo) = AX) + BU) 
AXbu)- AXhn)= Xo + BUlA) 


ALXu)-AXb)= 2 + BU) 


[nT-A]Xo)= X0 + Bula) 


| X (4) = [nI-A xo + [pa] BU(A) 


Xi e traje daria du X2 tapa au 
lindo, it odds imail mto 
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Vol = CXlol + DU(») 


Yol= €[nEAT'Y, + Chl-4] But) + DU(A) 


Vlo)= C [aI- AJ +[c [Zn] B +D] U(») 
ue pd) 


A 
Vo respolrcom  Yoç: reprta Com âm 
evdasda pala dica A'wtisk pula. . 
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O —. 


Yu)=|cCaT-A]'B + DJ]U(») 


for din bass, Seberm Ga 
To)= 6h) Vla) 
dum ds, 


Gln)= clar-a]"B+D 


CD —— 
O ——— 


ragga sileolo (e de pai ala ). 
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se falula à fracas du Lisusfiihnid 


Resolução da equação de estado no domínio temporal 


x = Ax+ Bu, x(0) = x0 


Parte homogénea n = Ax, x(0) = X, 
Exemplo 

x = ax, x(0) = x, (uma variável de estado) 
x(t) = e” x, de facto, 


x = x,ae” = a(x,e") = ax 
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...e se À for uma matriz ? 


x = Ax, x(0) = x, (n variáveis de estado) 
dD)=e"x, 2999 


e fé t* 
Pelada dd dd Bo 
21 3! k ! 


t? tê t* 
e =sASADrA RA dA Ea 
2 3! k ! 
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Matriz de transição de estado, D(t) 


x(t) = D(t)x, 


Conhecida a condição inicial, para calcular o estado em 
qualquer instante futuro basta multiplicar o estado inicial pela 
matriz de transição de estado nesse instante; isto é, esta 
matriz transita o estado inicial no instante O para o instante t, 
e daí o seu nome. 


No caso linear invariante (coeficientes constantes) 


XD) = ex, D(t) = e” 
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Polinómio característico da matriz A 
21-A/=0 
QM)=1"+a, A" +.+aA+ a =0 


Valores próprios de A: raízes da sua equação característica 


Vectores próprios à direita de A: vectores v tais que 
Av=AveS 4v- Av=0<>[41-Alv=0 
Vectores próprios à esquerda de A: vectores w tais que 


w1i=wAsSwiA-wA=-06w[41-4]=0 
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Estrutura própria da matriz A (caso de valores próprios distintos) 


Conjunto dos 
Valores próprios A, 


Vectores próprios 


Dr 


direita v,; 


a 


esquerda w, 


Matrizes 
Constituintes 
de A 


A matriz A pode ser reconstruída a partir da sua estrutura 
própria (caso n valores próprios distintos) 
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Pelo Teorema das Matrizes Constituintes 
- Att 
e” — ; ervw 
JJ 
e! 
e portanto 


x(t) = Ser v; W; |x(0) 


=evw x0)+et.y wx(0).+..+e.v w.x(0) 
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Este resultado diz-nos que, para uma dada condição inicial 


LÁ 


- a trajectória temporal do estado é uma soma 
ponderada de exponenciais dos valores próprios da 
matriz A ( os termos e). 


* São Os vectores próprios que estabelecem os 
coeficientes de ponderação. 


| a Ai 
Se existe um 4 real positivo lime” =o00 
ts 00 


...€ O sistema é instável em relação às condições iniciais!!! 
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Caso de um par de valores próprios 
complexos conjugados: 


e Se a parte real é negativa, contribui com uma 
sinusóide decrescente 


- Se a parte real é positiva, contribui com uma 
sinusóide crescente. 


- Se a parte real é nula, contribui com uma 
sinusóide de amplitude 


- constante se for simples, 


- crescente se for múltiplo. 
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Condição necessária e suficiente de 
estabilidade às condições iniciais: 


e Que todos os valores próprios de A 
tenham parte real negativa ou nula, 
mas neste caso têm que ser 
simples (multiplicidade 1). 
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Solução da equação de estado completa 


x = Ax+ Bu 


x(t) = e“ x(0) + [e Bu(o)d T 


O 


Qua nus 


X(D)=X,+X. 
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Pontos singulares de sistemas não-lineares 
Para um sistema não-linear definido por 
x= [(x,u) 


São pontos singulares (também chamados singularidades ou 
estados de equilíbrio) aqueles em que 


x=06 f(x, u)=0 
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Exemplo 


'=x1+x2-2u 
x2'=x2 (x1 - x2) 


x1 


traçado com o pplane.jar disponível em 


2/08/2022 


Www.cs.unm.edu/—joel/dfield 


https: 
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Linearização em torno dos estados de equilíbrio 
(também chamados singularidades) 


Série de Taylor 


[ (XxX + AX, + AM,...X o + AX ou, + Au,,U,, + Ay,,...,U, + AU o 


fi(x Xts> ATER -Anso ço 289º AD o 5 RA, 


+ termos de ordem superior 
Gi (x, + AX, >, + Ax,,...X + AX ou, + Au,,u,, + Au,,...,U, + Au )= 


= 09, — 09, 
S(Xiso X350 +: My5o Uso Laço ces Um ) FD AR, +> ay, 


k=1 Xi k=1 U, 


+ termos de ordem superior 
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Definindo os J acobianos 


Eae ho 
OX, OU, Ou, 
= B=|... 
EM fo df, 
OX, E OU, ou, E 
og 
ox, dg dg 
Ou, Ou, 
o D=) ... aid 
og, og, og, 
OX, a ou, Qu, a 
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Substituindo na série de Taylor 


Xs — f(xu,) y, = 9(Xu,) 
(xs + Ax) = f(x, + Ax,u, + Au) Y. + Ay=g(x + Ax,u, + Au) 
xs+Ax= [(x,u,)+ AAx+ BAu y.+Ay=g(x,u)+CAx+DAu 


fo = AA Dn 
Ay = CAx + DAu 
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Exemplo 1 


2 
à us o : u 
u,=1 
X, 22, + AX (11) O X, —2.X, + AX [2,0 


C=|x Xhn D=[0] C=|x Xloo D=[0] 
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C=[1 1] 


eig(A) 
ans = 


-1.4142 
1.4142 


D=[0] 
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2 
2 
XxX. = u 
S 0 S 
N geo do 
0 Io2 lo 
C=[2 1] D=[0] 
» eig(A) 
ans = 
À 
2 


áve 


t 


LÁ 


o Ins 


As singularidades são um ponto sela e um n 


x1+x2-2u 
x2 (x1 - x2) 


x1' 


x2'= 


-0.5 


traçado com o pplane.jar disponível em https://www.cs.unm.edu/-joel/dfield/ 
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Exemplo 2 


x =(x, ID) -—u y 


| 
bas 
IN9 


x =(%-2)-u 


(,1)-u=0 SM 1 =usS(x -D=+vu 
(4-22-u=0 S(y-2 =uS(x,-D=+Vu 


S (x, D=+Vu O x,- ENA 
S (-D=+Vu S x =2+u, 
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o 


nó instável 


o SJfao poi 4 


ponto sela | ponto sela nó estável 


(DADC/DEI/FCTUC/LEB/MCPF/DSBF/Cap. 4 219 


traçado com o pplane.jar disponível em https://www.cs.unm.edu/—joel/dfield/ 
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Conclusão 


A representação de estado aplica-se de igual modo aos 
sistemas lineares e não lineares, variantes ou invariantes. 


No caso linear obtém-se uma representação matricial. As 
propriedades dinâmicas do sistema são dependentes dos valores 
próprios da matriz de estado, tal como são dependentes dos polos 
da função de transferência na representação no domínio complexo. 


Os sistemas não lineares podem ter zero, um ou vários estados de 
equilíbrio, ou singularidades, ou pontos singulares, para a mesma 
entrada. Alcançam um ou outro conforme as condições iniciais. 


Aproximando as funções de estado e de saída pela série de Taylor 
nos pontos singulares, desprezando os termos de ordem superior 
a primeira, obtém-se um sistema linearizado. A estabilidade local, 
em torno da singularidade, é definida pela estabilidade do 
sistema linearizado. 


(DADC/DEI/FCTUC/LEB/MCPF/DSBF/Cap. 4 221 


Bibliografia 


Baura, G. D., System Theory and Practical Applications of Biomedical Signals, , 
John Wiley and Sons, 2002 


Bruce, Eugene N., Biomedical Signal Processing and Signal Modelling, John 
Wiley and Sons, 2001. 


D'Azzo, John J. and Constantine H. Houpis Linear Control System Analysis and 
Design with Matlab, 5t" Ed. Marcel Dekker 20083. 


Carvalho, M., Sistemas de Controlo Automático, LTC- Livros Técnicos e 
Científicos Editora, 2000. 


Chen, T. Systems and Signals Analysis, 2nd Ed, Sunders College Publ, 1994 


Hoppensteadt, F.C. and Charles S. Peskin,, Modeling and Simulation in 
Medicine and the Life Sciences , Springer Verlag, 2000 


Ribeiro, M. Isabel, Análise de Sistemas Lineares,., IST Press 2002 


pplane is copyrighted by 


john C. Polking, Department of Mathematics, Rice University. 
Mais informação em : https://www.cs.unm.edu/—joel/dfield/ 5/08/2022 


Versão atualizada para Matlab 


https://www.mathworks.com/matlabcentral/fileexchange/61636-pplane 5/08/2022 


(DADC/DEI/FCTUC/LEB/MCPF/DSBF/Cap. 4 Zoo 


Capítulo 5 


COMPORTAMENTO DOS 
SISTEMAS NÃO LINEARES 
E CAÓTICOS 


Pontos singulares (singularidades) de sistemas lineares 


Todas as curvas de fase deste capítulo foram traçadas com o pplane.jar disponível em https:/Awmww.cs.unm.edu/-joel/dfield 


Pont sela 


v'=-(kx+dv)im 


a E = 
ea x 
A e 


>» 0--— 


Foco estável (poço-sink, atractor) 


Nó estável (poço-sink, attractor): 


1 vp real positivo + 
1 vp real negativo 


2 vp complexos conjugados 


com parte real negativa 


2 val. Próp. (vp) reais negativos 


Foco instável (fonte-source, repelling): Um centro 


Nó instável (fonte-source, repelling) 


2 vp complexos conjugados 2 vp imaginários conjugados 


com parte real positiva 


2 vp reais positivos 
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Modelo de Lotka-Volterra 


04 b=0.01 


a 
(o) 


(a - b predador) presa 


predador ' 


presa 


0.005 
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Estável ou instável ? 


Estabilidade local 
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Incerteza nas condições iniciais ... 


... Impossível predizer o comportamento futuro 
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Edward Norton Lorenz em 1992: 
“Predictability: Does the Flap of a 


Butterfly's Wings in Brazil set off a 
Tornado in Texas ?” 
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Oscilador de van der Pol: modelizar a actividade 
oscilatória do coração humano 


db dy-Xs 
do |! 
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Sem excitação externa dx e 


seda oba alado ] dt Ê, 


“um ciclo limite estável (limit-cycle) 
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Com 


excitação externa 
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«. princípio do pacemaker 
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Sistemas oscilatórios não lineares 


se excitarmos o sistema com uma frequência muito afastada da sua 
frequência natural, o comportamento obtido resulta em geral num 
sinal aperiódico. 


se a frequência de excitação se aproxima da frequência natural, o 
oscilador “adota” a frequência de excitação. 


Esta propriedade é muito importante na prática: ela está na base dos 
pacemakers cardíacos. 


Também se verifica em : 

-adaptação dos ritmos biológicos ao ciclo dia-noite, 

-acoplamento entre a respiração e a pressão sanguínea, 

-Sincronização dos geradores fisiológicos de padrões centrais 
durante a caminhada e a corrida. 


Outros fenómenos fisiológicos podem ser modelizados com recurso a 
osciladores não lineares. Por exemplo o oscilador de Poincaré é 
apropriado para simular as arritmias cardíacas. 
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Populaçao 


Bifurcações e caos 


Uma espécie biológica com efeito limitante: 


Xe — Ax dx, )= Alx, = Ryo), xe[0,1] 
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A=3,5 


Populaçao 
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O O e O O 


dia 
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Bifurcações 
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O delta de Feigenbaum 


5 = Ada 
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Exemplos de sistemas fisiológicos com comportamento 
caótico 


Regulação da densidade de neutrófilos no sangue 


Variabilidade cardivascular 
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Conclusão 


Os sistemas não lineares permitem exprimir características 
importantes dos processos fisiológicos. 


O plano de fase (sistemas de 22º ordem) permite analisar a 
estabilidade de cada ponto singular. 


Uma bifurcação é um ponto (no espaço dos parâmetros ) em que 
o período de oscilação duplica. 


O comportamento caótico (caracterizado por ser aperiódico) 
está associado a sucessivas bifurcações, que definem o caminho 
para o caos produzido por variações de parâmetros do sistema. 


O número de Feigenbaum, considerado uma constante universal, 
permite encontrar as bifurcações. 
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Tudo o que estudámos é fundamental para modelizar 
processos fisiológicos: 


* Analogias e sistemas análogos 

e Obtenção das equações diferenciais 
* Função de transferência 

* Espaço de estados 

* Curvas de fase 


* Sistemas não lineares e comportamento 
caótico 


* Implementação em Matlab/Simulink 
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